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Abstract 

Assuming only the non-zero electron and tau neutrino components eee, Cerj Crr of the non- 
, standard matter effect and postulating the atmospheric neutrino constraint e^r = |eeTp/(l + fee); 

pi ^1 we study the sensitivity to the non-standard interaction in neutrino propagation of the T2KK 

^ \ neutrino long-baseline experiment. It is shown that T2KK can constrain the parameters |eee| ^ 1, 

leerl ~ 0.2. It is also shown that if legrl and ^13 are large, then T2KK can determine the Dirac 
phase and the phase of separately, due to the information at the two baselines. We also provide 
an argument that the components \eaii\ {a = e, /i, r) must be small for the disappearance oscillation 
r~ I ! probability to be consistent with high-energy atmospheric neutrino data, which justifies our premise 

Qh| that these quantities are negligible. 
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I. INTRODUCTION 



It has been shown by experiments with solar and atmospheric neutrinos yy] that neutri- 
nos have masses and mixings. In the standard three-flavor framework, neutrino oscillations 
are described by three mixing angles, 9i2, O13, 623, one CP phase 6, and two independent 
mass-squared differences, Am|^ and Amg^. The values of the set of the parameters (Amgi, 
612) and (Amgj^, ^23) have been determined to a certain precision by solar and atmospheric 
neutrino experiments, respectively. On the other hand, only the upper bound on ^13 is 
knowiJ3, sin^ ^13 < 0.04 at 90%CL, and there is an absence of information on 5. In future 
neutrino long-baseline experiments (see, e.g., Ref. [sl), the values of 6*13 and 6 are expected 
to be determined precisely. As in the case of B factories |9|,[l0|, such highly precise measure- 
ments will enable us to search for deviation from the standard three-flavor oscillations. One 
such possibility, which will be discussed in this paper, is the effective non-standard neutral 
current-neutrino interaction with matter (ll|-[l3| , given by 

(1) 



where / and /' stand for fermions (the only relevant ones are electrons, u, and d quarks), 
Gp is the Fermi coupling constant, and P stands for a projection operator that is either 
Pl = {I — 75)72 or Pji = {1 + 75)72. In the presence of the interaction Eq. ([1]), the standard 
matter effect Il4| is modified. Since we discuss the long-baseline experiments on the 
Earth, we make an approximation that the number densities of electrons (iVg), protons, 



and neutrons are equalJl By introducing the notation e^/? 
hermitian 3x3 matrix of the matter potential becomes 

f 1 + eee e, 

A = A 



eP 



+ 3e^J + 3e 



dP 



V 



the 



(2) 



where A = y/2GFNe. 

Constraints on e^^ from various neutrino experiments have been discussed in Refs. 15 
Since the coefficients e 



23 



in Eq. 



21) are given by e^/j ~ e^^ + 3eJ^^ + 3e„^, considering 



the constraints by Refs. 15-23|, we have the following constraints 24] at 90%Cljl: 



|e,J < 4 X 10° 



~efi I 



'-MM 



< 3 X 10- 

< 7 X 10 



^fir I 



(3) 



< 3 X 10° 

< 3 X 10-^ 

< 2 X 10^ 

From this prior study, it is known that the bounds on eee, Cer and e^r are much weaker than 
^afi (« = e, fi, r), and typically ~ (a, /3 = e, r) is allowed. 

On the other hand, the new physics with components (a, /3 = e, r) should be con- 
sistent with the high-energy atmospheric neutrino data, which suggest the behavior of the 



In Refs. , a global analysis of the neutrino oscillation data has been performed, in which a non- 
vanishing best-fit value for 613 is obtained. This result, however, is compatible with ^13 = at less than 
2cr, and it is not yet statistically significant enough to be taken seriously. 

This assumption is not valid in other environments, e.g., in the Sun. 

Here, we adopt the conservative bounds on eap which were derived without using the one- loop arguments, 
because the bounds obtained by the loop-contributions are known to be model-dependent [25[ . However, 
if we accept results based on one-loop argui nents, then we get the following bounds at 90%CL: 



P|2 



and |ef^| < 3.1 x 10 
\et^\ <0.32 



I3e: 



dP 



27 



27 



26|428|, where \el^\ < SxlO-^llSl, |e^^| < 3.1x10 



m, \e^P\ < 0.28 



27 



28| and |eff | < 0.28 
2 



27, 



27, 



11/2 ^ Y in-3 

1 are used; {e^rl < Ep(kef Jl+l3e^r^|' + - 1.7 |27|, l28|, where 



281 are used. 



disappearance oscillation probability 



sin^ 26 



atm 



sm 



V ^ 



oc 



1 



(4) 



atm and Am^^jjj are the oscillation parameters in the two-flavor formalism. Note 



where sin^ 26, 

that the terms of 0{E'^) and 0{E~^) are absent in Eq. (jl]). As shown later, the elements 
Cq/x ~ ''") should be small so as to produce no term of 0{E^) in Eq. (jl]). Furthermore, 



absence of terms of 0{E ^) in Eq. (j4]) implies 



:i + ee 



(5) 



as pointed out in Ref. 
matrix (12]) with e 







30j . When Eq. is satisfied, two of the three eigenvalues of the 
e, fj,, t) become zero. Only in this case, one of the frequencies of 



the three oscillation modes at high energy coincides with the one in the vacuum oscillation, 
and the disappearance oscillation probability of z/^ behaves as in Eq. (jl]). The effect of 
the non-standard interaction in propagation for solar neutrinos has also been discussed in 



and [32[ give a constraint —0.06 < e'^^ sin ^23 < 0.41 



Refs. [19|, [221, |23|, jSl], |32l , and Refs. |2 
(at 90%CL) and |ef]f | < 0.4 (at Ax^ = 4 for 2 d.o.f.), respectively. 

The sensitivity of the ongoing long-baseline experiments to the non-standard interaction 
in propa gati on has been studied for MINOS 33 in Refs. [33^37j. and for OPERA [38*] in 



On the other hand, for the sensitivity of future long-baseline experiments, 

an extension of the T2K 



Refs. 

Ref. [41[ provided the sensitivity of the T2KK experiment |42l, 143 

neutrino oscillation experiment [3] with a far detector in Korea, in the two-flavor analysis 
with = {a = e, fi,T). The sensitivity of the reactor and super-beam experiments was 
discussed in Ref. (isl . and the sensitivity of neutrino factories js], |46| has been discussed by 
many authors |47H54l| . 

In the present paper we analyze the sensitivity to the parameters e^/s of the T2KK 



experiment, assummg e 







a 



e, At, r) and err = l^eTp/ (1 + Cee)- We do not discuss the 



so-called parameter degeneracy [55H58| , since little is known about parameter degeneracy in 
the presence of the new physics, and the study of the subject is beyond the scope of this 
paper. 

The paper is organized as follows. In sectHIl we discuss the constraints on the 
parameters from the high-energy atmospheric neutrino data. In sect lllll we analyze the 
sensitivity of T2KK to the new physics parameters. In sect Jl VI we draw our conclusions. 
In the appendices RHHl we provide details of the derivation of the analytic formulae for the 
oscillation probabilities and their high-energy behaviors. 



II. ATMOSPHERIC NEUTRINOS AND THE NON-STANDARD INTERACTION 
IN PROPAGATION 



In Refs. [29], |30| , it was pointed out that the new physics should be consistent with the 
constraint imposed by the atmospheric neutrino data, which suggest that the disappearance 
oscillation probability at high-energy satisfies Eq. (jlj. In the standard three-flavor scheme, 
the high-energy behavior is 



~ lim 



m3 I 



1 



m3 I 
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sin^ 2^23 
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|2„2 



sm 



/AL 



S23 sin^ 2^13 sin' 



/AL 



(6) 



where AEj^ = Ej — i^^ ^ Am'jf^/2E. (See the appendix |X] for derivation of the oscillation 
probability in constant density matter and the appendix [B] for derivation of its high-energy 
behavior (jS])). In the following discussions, we assume the generic matter potential ([2]), and 
derive the high-energy behavior 



1 

^2 



(7) 



and postulate |co| ^ 1 and |ci| <C 1. Note that the term of corresponds to the 

standard results (jl]) in the two-flavor scheme or (jH]) in the three-flavor case, where the 
information on the atmospheric neutrino oscillation parameters appears. In the presence of 
the matter potential ([2]), as we discussed in detail in the appendix O, assuming the constant 
density of matter, we can obtain the analytic form for the behavior (j4]) of the disappearance 
oscillation probability at high energy, in which \A/AE3i\ = V2GFNj{\Aml-^\/2E) < 1 is 
satisfied. The leading term cq in Eq. ([7]) is given by 



Co 



1 + eee + ej 
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^1 ~l~ ^ee) l^/xr 



X sm 
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2Re(6e^e^T- (Cer J 



1 + e.e + 
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-|- 67-7- 1 Cg^ I 2Re(6e^€^T- (€er) ) 



+4 



-e/i I 



+ e 



MM I 



+ e 



MT| 



1 -1~ 6gg ~1~ €7-7- 

(1 + eee)|e^^| 



+ e 



rr I ^e/i | 



^1 -|- 6ee ~l~ Crr) 
/ {1 + e^e + err)AV 



1 + eee + e-r 



X sm 



(8) 

(see appendix 



For Eq. ([8]) to be consistent with Eq. (jl]), we should have 
p. 

In Refs. [igI . 18, 59|, the two- flavor analysis of the atmospheric neutrino data with the 
matter effect 



A 



was performed. In this case, in the limit E ^ 00 the disappearance oscillation probability 
results in 



1 - P(z/^ u^) ~ 







2 


(c/./. - e,,)V4 + 




2 



sin' ALJie 



"-MM 



This suggests that the two parameters \e^r\ and |e^^ — e^rl should be small so as to be 
consistent with the high-energy behavior (jl]), and in fact, the authors of Ref. 59|] obtained 
the bounds |e^^| < 1.5 x 10~' and |e^^ — e^^| < 4.9 x 10"^ at 90%CL. In the three flavor 
framework, as seen below, ea/s ~ C'(l) («,/? = e, r) can be consistent with the atmospheric 



4 



neutrino data as long as e 



3^(1 + 



is satisfied 29, 30 . We see that the two 



flavor constraints in Refs. 16|, uM, l59| are consistent with three flavor ones as follows: As for 

• I <^ 1 in the two flavor framework implies | e^r | ^ 1 • 



smce e 



<^ 1, the bound |e^^ 



The two flavor framework can be regarded as a subset of the three flavor case in the limiting 
case = ^er = ^13 = ^^^21 = O5 SO the constraint Eq. (E]) in the two flavor case leads to 
|e^T-| ~ 0. On the other hand, the bound on |e^^| in the three flavor case is independent of 
other components eai3, so the bound |e^T-| < C(10~^) is expected to be valid both in the two 
and three flavor cases. While |e^^| -C 1 and je^,-! ^ 1 were flrst shown in Refs. |l5| and 



16|, 

respectively, if we do not accept the one- loop arguments [l5j to constrain ee/j. as in Ref. 25[, 
then the observation in this paper that |ee^| ^ 1 follows from the the atmospheric neutrino 
constraint is new, although our discussions are based only on an analytical treatment El 

0. Then 



In the following discussions, therefore, let us assume that e 



we obtain the coefficient Ci of the term of the next-to-leading order in 1/E in Eq. ([7]) (see 
appendix [D|) : 



Cl 



(1 + e,e)AAts2 




(9) 



where 



At32 = 



2(1 + eee) 
1 

2(1 + eee) 



1 + eee + err + J (l + eee + e,,)^ - 4(1 + eee)2C' 



;i + eee + e,,)2-4(l + eee)V- 



Again, for Eq. (jH]) with Eq. to be consistent with the high-energy behavior (jl]), we should 
have 



eer 



^TT , 

1 + eee 



0. 



This agrees with the conclusion e^r — leerP/ (1 + eee) in Refs. [29|, |30 



Thus, let us assume err — |eerP/(l + eee) = 0. In this case, it is convenient to introduce 
the new variable 



tp = tan (3 = 



1 + ee 



Then, we have the following high-energy behavior (see appendix lEl): 



(10) 



(AE3lM)2 



4. ^23 



1 - 



sin^ 20" 




(11) 



251 appeared, the bounds |e{^| < 0(10"'') (90%) in Ref. [151 based on the one-loop arguments 



Before Ref. 

were widely accepted, and this was used to justify the assumption e^^ ~ in Refs. 



29 



3Q|. 



5 



where c'(^ = cos9'(^, Cfs = cos/3, Sjs = sin/3, and the new angle 9'{^, which is introduced in the 
appendix [F] to diagonahze the mass matrix in the presence of the new physics, is defined by 
Eq.(lF6]). 

Comparing Eqs. and (fTTj) . we see that if 

1 + fee = 4 (12) 

is satisfied, then by introducing the two effective mixing angle 

sin^^-^^ (13) 

^^3 = O'l, , (14) 

Eq. (fTT]) shows almost the identical behavior as that of the standard scheme ([6]). 
A few remarks are in order. 

Firstly, although Eq. f ll2p is satisfied only in a narrow regionjl as long as Eq. (1121) holds, 
the high-energy behavior of the disappearance oscillation probability coincides with that 
of the standard three-fiavor scheme. Off this upper half circle, equivalence between the 
behaviors of Eqs. f|TT]) and ([6]) are lost, but it is expected that due to the experimental errors 
around this upper half circle, there exist some regions in which the behaviors of Eq. f|TT]) and 
([6]) are similar. 

Secondly, Eq. f|T^ indicates that the angle 6'(^ plays a role similar to that of ^13 in the 
standard scheme. Note that the corrections in Eq. (jH]) due to 6*13 were not discussed in 
Refs. [2^, lO] , where it was suggested that the quantities that appear in Eq. flTT]) imply the 



effective two-flavor mixing angle, sin 6'atm = (1 + ^^)'523/(l + ^23^) ^^"^ effective mass- 
squared difference Aml^^^ = Am|2(l + -523^1) /(I + While the former is exactly the same 
as Eq. (fT3|) in the limit 6*13 — )■ 0, the latter does not appear in our result. This is because 
the correction factor (1 + 523^^) /(I + 1|) naturally arises from the three-flavor contributions, 
i.e., from the 613 dependent terms, and we need not normalize Am^. If we postulate the 
effective mixing angle to be ^23 = in Eq. f|T3l) . then C23 = cos 6^23 > can be expressed 
by /3 and 613 as 



SpCpSisCisCOS^ ^ ] l + clsl^ ^ / g/jC^gl3Ci3COS'l> \% ^^^^ 
2 — S^Ci3 i 2 ~ '^§'^13 V 2 — ■5^Ci3 

where we have introduced 

$ = 5 + arg(ee^). 

In the limit 613 — )• 0, Eq. (ITHj) agrees with the expression C23 = 1/(1 + c^) obtained in 
Refs. MS. 



1/2 



Thirdly, Ref. [GOj performed a three-flavor analysis of the atmospheric neutrino data, and 
the authors concluded that the atmospheric neutrino data alone gives sf^ < 0.14 (0.27) at 
90%CL for a normal (inverted) mass hierarchy. This implies that the range 5^3 < 0.14 is 
consistent at 90%CL with the high-energy atmospheric neutrino data, i.e., the upward going 
H events. In the present case, we found that the value of {si^)"^ can be made smaller than 



^ In fact, Eq. (fT2|) stands for an upper half circle around a center (-1/2,0) with a radius 1/2 in the (tee, Icerl 
plane. 
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0.14 in almost all the region for < sin^ 2^13 < 0.15 and < s| < 0.5, by adjusting the value 
of $. With the conditions (fT2|) . (fT3|) . ([HI), therefore, the region around the upper half circle 
for |eer| ~ 0.5 and —1/2 < 1 + eee ~ is expected to be consistent with the atmospheric 
neutrino data. 

Thus, taking into account the various constraints described above, we will work with the 
ansatz 



/1 + 



A = A 













V(l + eee); 



(16) 



in the following discussions. 



III. SENSITIVITY OF T2KK TO e^e AND eer 

In this section we discuss the sensitivity of the T2KK experiment to the non-standard 
interaction in propagation with the ansatz f lT6|) . Since err is expressed in terms of e^r and 
eee, the only new degrees of freedom are eee, |eer| and arg(eeT)- Firstly, in sect. lIII A| we 
briefly describe the setup of the T2KK experiment. Secondly, in sect. lIII B| we consider the 
(eee, kerl) plaue and discuss the region in which T2KK can discriminate the non-standard 
interaction in propagation from the standard three-flavor scenario. Thirdly, in sect. lIII C| we 
study the case in which new physics can be discriminated and discuss how precisely T2KK 
can determine eee and |eer|- Then, in sect. lIII D| we consider whether the two complex phases 
5 and arg(eer) can be determined separately. 



A. The T2KK experiment 



The T2KK experiment |42l . |43| is a proposal for the future extension of the T2K experi- 



ment j4J]. In this proposal, a water Cherenkov detector is placed not only in Kamioka (at 
a baseline length L = 295 km) but also in Korea (at L = 1050 km), whereas the power of 
the beam at J-PARC in Tokai Village is upgraded to 4 MW. As in the T2K experiment, it 
is assumed that T2KK uses an off-axis beam with a 2.5° angle between the directions of the 
charged pions and neutrinos, and the neutrino energy spectrum has a peak approximately 
at 0.7 GeV. Because the two detectors are assumed to be identical, some of the systematic 
errors cancel. Also, because the distances of the two detectors from the source are differ- 



ent, parameter degeneracy in the three-flavor oscillation scenario |55l-l58| is expected to be 
resolved with this complex 42|, |43|, I6ll-l64 , 

In our analysis, we 



In this paper, we assume the same setup as that in Refs. j41|, |62 



use the disappearance channel z/^ — )■ i/^ and z/^ — )• z/^, the appearance one — )■ Ue and 
i^fi i^e, and data from single- Cherenkov- ring electron and muon events. We assume that 
the measurement will run for 8 years in total, 4 years each for the neutrino and anti- 
neutrino beams. The fiducial volume of each detector is 0.27 Mton. The density of the 
Earth is assumed to be p = 2.3 [g/cm'^] in the case of Tokai-Kamioka, and p = 2.8 [g/cm^] 
in the case of Tokai-Korea. The electron fraction is assumed to be 0.5. The energy 
resolution is considered to be 80 MeV. We use various information such as the neutrino flux 
from Ref. 65 . 
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B. Bounds 

Firstly, we discuss the case of the region (eee, l^eA), in which we can distinguish the new 
physics with ansatz flTBl) from the standard three-flavor scheme. To perform such a test, we 
introduce the following quantity: 



Ay^ = min 

par am, r 



' ^ 1-^=1 i=l '-'i \^) 1=3,7 



'=1,2,7 



20 



j=l '^i 1=4,5,7 



1=4,6,7 



+ E(-)' + AX^nor 
Z=l ^« 



(17) 



where the prior Ax^^ioj. is given by 



Ax 



prior 



2.7 X (sm^2e^ 



{omy 



23 -1.0)^ ^ (sin^gi3-0.02)^ 



(lAmjil -2.4 X 10-=^[eV'])^ 
(1.5 X 10-4[ev2])2 



(0.01)^ 



In principle we could perform an analysis without the prior AXpj.iQj,, but in that case it 
would take more computation time by minimizing Ax^ for the parameter region which 
is already excluded by the present data of the atmospheric and reactor experiments. So 
we have included the prior in our analysis to save computation time. In Eq. f|T7|) . N^{e), 
A^°(/i) (i?°(e), and B^{fi)) are the expected signal (background) numbers of events in the 
presence of the new physics f|T6|) . while iVj(e), Ni{fi) {Bi{e), and Bi{fi)) are the expected 
signal (background) numbers of events in the three-flavor framework with the standard 
matter effect. All these numbers except Bf{e) and -Bj(e) depend on the neutrino oscillation 
parameters. The indices i and = 1, ■ ■ ■ , 4 stand for the number of the neutrino energy bin 
for electrons and muons and the four combinations of detectors in Kamioka and Korea with 
the neutrino and anti-neutrino beams, respectively. For the electron events, there are five 
energy bins (400-500 MeV, 500-600 MeV, 600-700 MeV, 700-800 MeV, and 800-1200 MeV), 
whereas for the muon events, there are twenty bins from 200 to 1200 MeV with 50 MeV 
width. o"j(£) {i = e, /i) stands for the statistical uncertainties, whereas ee (£ = 1, ■ ■ ■ , 7) 
stands for the systematic uncertainties in the expected number of signals and backgrounds. 
Ax^ is defined by minimizing the quantity inside the square bracket in Eq. (1171) with respect 
to the uncertainties as well as the oscillation parameters (|Am|]^|, sign(Am|]^), 6*23, S) of 
the standard three-flavor scheme, on which the numbers of events Ni{e), Ni{^), Bi{e), and 
Bi{ii) depend. The uncertainties in -Bj(e) and A^i(e) are represented by 4 parameters ej (j = 
1, 2, 3, 7). The backgrounds in the muon events are referred to as non-quasi-elastic events in 
The uncertainties in Bi{fi) and Ni{fi) are represented by 4 parameters ej [j = 
7). The parameter /(e or /i)* indicates the possible dependence of the parameter ej on 
the i-th energy bin. ei stands for the uncertainty in the overall background normalization for 
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electron events with (Xi = 0.05. €2 is the energy- dependent uncertainty for the background 
electron-hke events with a function /(e)^ = {{E^ - 800 MeV)/400 MeV) and cts = 0.05. €3 
is the uncertainty in the detection efficiency for the electron signal events with = 0.05. €4 
is the energy- dependent uncertainty for both the muon signal and background events with 
the function /(/i)^ = {E^ - 800 MeV)/800 MeV and = 0.05. €5 is the uncertainty in the 
signal detection efficiency for the muon signal events with = 0.05. eg is the uncertainty 
in the separation of quasi-elastic and non-quasi-elastic interactions in the muon events and 
o"6 = 0.20. ey stands for the uncertainty in the neutrino flux in Korea, and 5"7 is assumed to 
be the predicted flux difference between those in Kamioka and in Korea, given in Ref. (66| . 

In Eq. ( !T7|) . the numbers of events iV°(e), iV°(/i) (i?°(e), and _B°(/i)) depend not only on 
the new physics parameters eee, l^eA, arg(eer) but also on the standard oscillation parameters, 
which we denote as 9i2, ^13, ^23, ^^21^ ^^315 ^- Here, we take the best-fit values for 
most of the standard oscillation parameters as the reference values: 

sin2(2^i2) = 0.87 

sin2(2^23) = 1.0 

Am^i = 7.9 X lO-^eV^ 

Am^2 = 2.4 X lO-^eV^ (18) 

On the other hand, since we have no information on ^13 and 6, we will take several reference 
values for these parameters. 

The results are shown in FigsHHU where the curves are drawn at 90%CL (Ax^ = 4.6 
for 2 degrees of freedom). The new physics with the ansatz (fT6|) can be distinguished from 
the standard three-flavor scheme outside the curves. Four different choices for the phases S 
and arg(eeT) are taken, where the sum of the two phases is the same in each figure. It has 



been known |48l . |67| that in the limit of Am^i — ?• 0, the oscillation probability P{i^^ — )■ i^e) 
depends only on the combination 6 + arg(eer) of the phases. If the four curves in FigsilHl] 
coincided with each other, then it would mean that the contribution of the solar mass- 
squared difference were small. From these figures, however, the behaviors of the four curves 
are different even if (5-|-arg(eer) = constant, so the contribution of Amli is not negligible. This 
is because we are considering the oscillation probabilities in Korea, where Aml^L/4E ~ 0.3. 
The analytic form of the oscillation probability Pii^fi — ?■ i^e) is given in the appendix[Fl where 
the correction to ^(t'^ — i^e) to the first order in Amli is also given. The corrections to the 
energy eigenvalues Ej {j = 1,2,3) are proportional to either sin^ ^'/j or cos^6''/2, where 6'(2 
is defined in Eq. flF7p . and depends mainly on arg(eeT)- From these figures we conclude 
that the T2KK experiment can distinguish the new physics with ansatz f lT6|) at 90%CL 
approximately for \eee\ ~ 1 and |eer| ~ 0.2. In other words, if T2KK lacks evidence of a new 
physics, then with the ansatz f|T6|) . we can put bounds on the two parameters: \eee\ ~ 1 and 
leerl ~ 0.2. While the bound on eee is modest, the one on |eer| is impressive compared with 
the present bound (E]). 

C. Precision in determination of eee, Icerl 

Let us now turn to the case with affirmative results in the T2KK experiment, i.e., we will 
discuss the points that lie outside the curves in the (eee, l^erl) plane in FigsJTI- HlFI 



^ Since we have the fixed value of 6*13 , we do not have the 613 — Cer confusion in Ref. [49| , in which it was 
shown that the same neutrino survival probability is produced by the different pairs of (^13, Cer)- 



9 



0.25 



0.2 



arg(e )=0 5=0 



arg(E..)=lo 5=270 --- 

arg(e„>180 5=180 

arg(EeJ=270 5=90 — ■ 



0.15 



0.05 




0.25 



0.15 




0.05 



0.25 



0.2 



arg(eg,)=0 5=180 

arg(e„)=90 5=90 --- 
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FIG. 1. Region in which the new physics is discriminated at 90%CL from the standard three-flavor 

^13 



scheme for sin^ = 10 ^. 



First, we discuss the experimental errors in eee and |eer|, i-e., the correlation of errors for 
these two variables. In this case, we introduce the following quantity, which is similar to 
Eq.([I7D: 

AX^ = . min [ E { E ^T-dNKe) + B^{e) - iV,(e) ^ (1 + /(e)!^) 

' ' ^k=l ^ i=l "i 1=3,7 

-B^{e) E (l + /(e)jQ)r 

1=1,2,7 

20 1 

i=l \h^) 1=4,5,7 



par am, ( 



1=1 



(=4,6,7 



(19) 



where most of the definitions are the same as those in Eq. (fTTl) . The only differences between 
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FIG. 2. Region in which the new physics is discriminated at 90%CL from the standard three-flavor 
scheme for sin^ 2^13 = 10~^. 



f fT7|) and f fT9|) are that the prior is absent in the latter, and that Ni{t) {i = e, fi) and 
Bi{£) {£ = e,/i) in the latter (in the former) are the expected numbers of events in the 
presence of new physics with ansatz ( I16p (in the standard scheme), respectively. Namely, 
both N^{i) {£ = e,/i) and {B^{i) {i = e,fi) and are the expected number of events in 
the presence of a new physics with the ansatz ( IT6|) in Eq. (fT9|) . The number of events 
^ii^) = ^5/^) depends on the parameters of the new physics (eee, Icerl, and arg(eer)) as 
well as the oscillation parameters of the standard scheme 612, O13, O23, Am^i, Am|^, and 
5. We fix 

^ee; l^erl at somc poiuts outsidc the curves in FigsJlVHl and evaluate as a 
function of eee and leerl, which appear in the argument of Ni{i) {i = e, /z). For simplicity, 
we assume the central values given in Eqs. (ITK]) for 9i2, ^23, A'^ii; ^^li- For ^13, we 
take a few representative values sin^ 2^13 = 10~^, 10~^, 0.12. We assume normal hierarchy 
and fix the value of the phases as 5 = vr and arg(eer) = tt for simplicity. As for the variables 
in Ni{i) {i = e, /i), for simphcity we equate the variables 6*12, 6*13, 6*23, Am2i, '^ml^, S, and 
arg(eer) to Ou, ^13, ^23, Amji, ^^31, ^, and arg(eer) in N^{i) {i = e, /x), respectively. In 
this analysis, we do not introduce any prior because it will be difficult to estimate it in the 
presence of the new physics. Thus, in Eq. ( 1T9|) . we only minimize the quantity in the square 
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FIG. 3. Region in which the new physics is discriminated at 90%CL from the standard three-flavor 
scheme for sin^ 2^13 = 10~^. 



bracket with respect to the parameters = 1, ■ ■ ■ , 7), and we evaluate Ax^ as a function 
of the variables eee and |eer|- 

The results are shown in FigJSl where the contours around the true points are depicted at 
90%CL (Ax^ = 4.6 for 2 degrees of freedom) for three different values of ^13. The straight 
lines Iter I = 1-1 X |1 + eee| in Fig|5] stand for the approximate bound from the atmospheric 
neutrinos, which lead to cos^ /3 > 0.45 or | tan /3| ^ 1.1 [soj, and we have examined only the 
points below these straight lines. As seen from FigJSl the error in leerl is relatively smaller 
compared with that in e^e for all the values of ^13. The experimental error in increases 
for smaller values of |eer| and 6*13. For |eer| ^ 0.3, therefore, the possibility of eee = cannot 
be excluded by the experiment, while |eer| = can be for |eer| <L 0.2. For 1 + eee < 0, the 
correlation in eee and |eer| increases, and in this case, even if T2KK can discriminate the 
new physics from the standard scheme, the determination of these parameters is difficult. 
We have also analyzed other combinations of the phases 5 and arg(eeT) and they share the 
same features as in FigjSl 



12 



0.25 



0.2 



arg(e J=0 5=0 
arg(E..)=§0 5=270 
arg(e„>180 5=180 
arg(EeJ=270 5=90 



0.15 



0.1 



0.05 




-1.5 -1 -0.5 0.5 1 1.5 



0.25 



0.15 




0.05 



0.25 



0.2 



0.15 



0.1 



0.05 



arg(eg,)=0 5=180 
arg(e.J=90 5=90 
arg(e.,)=180 5=0 
arg(eex)=270 5=270 



0.25 




0.15 



-1.5 



-0.5 



0.5 



1.5 




0.05 



FIG. 4. Region in which the new physics is discriminated at 90%CL from the standard three-flavor 
scheme for sin^ 2^13 = 0.12. 



D. CP violating phases 

In the ansatz (fT6|) . there are two phases S and arg(eeT), and if a new physics exists at all, 
then it becomes important whether we can determine these two phases separately. Thus, we 
address this question next. In this case, we can use the same Ax^ in Eq. ( IT^ . but there are 
two differences between this and sect JIII CI Firstly, here we vary the variables egg and |eer| 
in Ni{i) {i = e,fi) and minimize the quantity in the square bracket in Eq. f|T9|) with respect 
to these two parameters as well as the parameters eg. Secondly, Ax^ is plotted as a function 
of the two variables 6 and arg(eeT) in Ni{i) {i = e, /x) here, while it is plotted as a function 
of the two variables eee and |eer| in Ni{i) {i = e,fi) in sect JIII Cl 

The results at 90%CL are shown in Figs. E] and [7] for (eee, |eer|) = (0.8, 0.2) and (2.0, 2.0), 
respectively. Since we are discussing the cases that can be distinguished from the standard 
scheme, relatively large values of |eer| are chosen in both examples. Again 6 = arg(eer) = tt 
is assumed for the true values of the phases. To clarify the roles of the detectors at the two 
baseline lengths, separate contours are given for the result from the detector in Kamioka, 
for that from the detector in Korea, and for that from the combination of the two. As 



13 




FIG. 5. Correlation in (cee, |eer|): The contours around the true points are depicted at 90%CL. 
5 = aig^ier) = vr is assumed. 



in the standard three-flavor case, if ^13 is very small, neither the detector in Kamioka nor 
the one in Korea can provide any information on 6. As the value of ^13 increases, the 
sensitivity to S of the detectors in Kamioka and Korea increases. For larger values of ^13, 
the sensitivity to arg(eer) depends on the value of |eer|- For larger (smaller) values of |eer|, 
sensitivity to arg(eer) is good (poor). These features can be understood qualitatively by 
looking at the T violating term in the oscillation probability -P(t'^ — > z/g) (see appendix 
Ofor details). In the present case, there are two sources for T violation, the standard one 
Im{y3'^(y2^'^)*}std and the extra one Im{y3^'^(y2'^*^)*}NP because of the new physics. From the 
explicit forms (1G3P and (lG4p . the ratio of the two terms is roughly given by (std)/(NP)~ 
\AE3i\\Ue3\/{A\eer\\Ur3\) ~ 10si3/|eer|. This implies that if sin^ 2^13 > (<) O{10-^) and 
if leerl = 0.2, then the contribution of 613 is large (small). If {eerl is very large, then the 
oscillation probability and the number of events becomes so large that the sensitivity of both 
the detectors to the two phases increases, as shown in Fig. [3 
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FIG. 6. Correlation in {5, arg(eer)) for eee = 0.8, |eer| = 0.2. The contours around the true points 
are depicted at 90%CL. 
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FIG. 7. Correlation in {S, arg(eer)) for fee = 2.0, |eer| = 2.0. The contours around the true points 
are depicted at 90%CL. 

IV. CONCLUSIONS 



In this paper we have studied the sensitivity of the T2KK experiment to the non-standard 
interaction in propagation with the ansatz f lT6|) . 



To justify our premise, we have provided an argument that e 



fia 



(a 



e, /i, r) must be 

small for the behavior of the disappearance probability to be consistent with the high-energy 
atmospheric neutrino data. Using the analytical form of the disappearance probability, we 
showed that e^Q, (a = e,yU, r) ^ as well as err — |eeTp/(l + fee) are necessary for the 
disappearance probability to be consistent with the high-energy behavior 1 — — )■ z/^) oc 
in the high-energy atmospheric neutrino data. This speculation should be verified 
explicitly by numerical calculations in the future. 

With the ansatz (fT6|l . we have looked for the region in the (eee, l^eA) plane in which 
T2KK can distinguish the new physics from the standard three-flavor scheme. At 90%CL 
T2KK can discriminate the new physics from the standard case for approximately |eee| ^ 1 
and |eer| ^ 0.2. These values can be interpreted as bounds for these parameters, if T2KK 
has negative evidence for a new physics. While the bound on eee by T2KK is modest, the 



bound on |eer| by T2KK is much stronger than the present one jeerl ^ 3, and the latter 
is also stronger than those by other on-going long-baseline experiments such as MINOS or 
OPERA: If sin2 2ei3 ~ 0.07, then MINOS will give a bound |eer| ^ 1 at 90%CL [H-iJ, 
and the combination of MINOS, OPERA and Double-CHOOZ gives a bound -2.2 (-2.5) 



< 
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eer ^ 0.6(0.5) at 95%CL for sin^ 2^13 = 0.05 (0.1)|39|. 

On the other hand, if the true values of the new physics parameters he much outside 
these bounds, then T2KK can determine the values of eee, Icerl, and arg(eer)- In particular, 
if the values of 6'i3 and |eer| are relatively large (sin^ 2^13 > 0(0.01), |eer| ^ 0.2), then 
we can determine the two phases 5, arg(eeT) separately. This determination is possible, 
because the oscillation probability in Korea with the baseline L = 1050 km receives a non- 
negligible contribution from the solar mass-squared difference Am'^i, and it induces terms 
that approximately depend only on arg(eeT)- 

Since our analysis is based on the ansatz ( IT6|) . all the results derived in this paper are 
approximate ones. Nethertheless we believe that our conclusions are valid. 

Long-baseline neutrino experiments with longer baselines (L ^ 1000 km) are sensitive to 
the matter effect because of the matter effect contribution appears in the form of AL/2 ~ 
L/4000 km in the argument of a sine function in the oscillation probability. They are, 
therefore, also sensitive to a non-standard interaction in propagation. Among long-baseline 
experiments with longer baselines, T2KK, for which the matter effect and the contribution 
from the solar mass-squared difference are smaller than the one from the atmospheric mass- 
squared difference (|A£'2i| ~ A <C |A£'3i|), is different from neutrino factories, for which the 
contribution from the solar mass-squared difference is smaller than the matter effect and the 
one from the atmospheric mass-squared difference (|A_E'2i| ^ A ~ |A_E3i|). Consequently, 
the features of T2KK can complement those of neutrino factories. While T2KK is known to 
be powerful in resolving parameter degeneracy in the standard three-flavor scenario, further 
studies on the new physics potential of T2KK should be pursued. 



Appendix A: Analytic treatment of the oscillation probability in constant density 
matter 

Throughout this paper, we assume that the density of matter is constant. In this appendix 
we derive analytically the neutrino oscillation probability in constant density matter. Let 
us start with the Dirac equation 

= (usu-^ + A)^, (A1) 

(Jjv 

where \E' = {u^,, z/^, u^-)'^ is the flavor eigenstate, 

S = dmg{E,,E2,E3), (A2) 

is a diagonal matrix with the energy eigenvalues in vacuum, and we assume the nonstandard 
matter potential A defined in Eq. (|2]). In practical calculations, it is easier to subtract the 
mass matrix USU~^ + ^ by which only affects the phase of the oscillation amplitude. 
In the following, therefore, by Eq. flA2p we mean 

S = diag {EuE2, E3) -Eil = diag (0, AE21, AE31) , (A3) 

where AEjk = Ej - Ek ^ Am|^/2E. 

The 3x3 matrix on the right hand side of the Dirac equation 
diagonalized as: 

U£U-^ + A=U£U~\ (A4) 



can be formally 



16 



where 



is a diagonal matrix with the energy eigenvalues in matter. 

As in the case of the oscillation probability in vacuum, Eq. (lAip can be integrated as 



Thus the oscillation probability P(z/q, — )■ z/^) is given by 



13a 



lfexp{-t£L) 
5«,-4ERe(xfxf*)sin ^ 



2 / ^EjkL 



j<k 



-2 ^ Im (x/'^Xf *) sin [AEj^L 

j<k 



where we have defined 



yafi _ TT TT* 



(A5) 



(A6) 



'jk 



E, 



Ekj 



and a, /3 = (e, /i, r) and j, = (1, 2, 3) stand for the indices of the fiavor and mass eigenstates, 
respectively. Eq. ( ]A5I) shows that we can obtain the oscillation probability if we know the 
energy eigenvalues Ej and X°^. The former can be obtained by the eigenvalue equation 
\U£U-^ + A-tl\ = 0, while the latter can be obtained by the formalism of Kimura, 



Takamura and Yokomakura 68, 69 



Let us briefiy see how we get X"^ |67l-l71| . From the unitarity condition of the matrix U, 



we have 



Jal3 



uu 



-1 



«/3 



(A7) 



Note that the quantity X"^ was defined in Eq. (]A6|) . Furthermore we take the {a,/3) com- 
ponent of the both hand sides in Eq. ( 1A4|) and its square: 



a/5 



aP 











aP 





a/3 



j:Ua,E]u;^ = T.E]x: 



2 -^0/3 
j ■ 



From Eqs. ( ]A7p and (lASp . we have 



/ 1 1 1 

El E2 E-i 
V E! El El 



al3 



(XI 

xf 



yap 

-'2 

yap 



where we have introduced notations {Yi^ = dap)'- 



a/3 
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aP 



(A8) 



(A9) 



(AlO) 



Eq. f lAgp can be easily solved by inverting the Vandermonde matrix: 



Xf 

xf 



1 



Ai?2lA-E32 



\ A£;3iAE: 



{E2ES 
(E1E2 



-{E2 + E,), 
-(^3 + ^1), 
-(^1+^2), 



1) 
1) 



\y 



I 

Q/3 



(All) 



^32 



Note that Y^^ = 1, 2, 3) in Eq. (lAlOp can be expressed by the known quantities, so that the 
oscillation probability f lASP can be expressed analytically in terms of the known quantities. 

To summarize this appendix, the following is a set of the procedures to obtain the analytic 
form of the oscillation probability in the presence of constant generic matter potential: 



(i) Obtain the roots Ej of the eigenvalue equation \U£U ^ + A — tl\ =0. 



ii) Obtain the coefficients Xf in Eq. fIXTT]) by evaluating Y^" = {USU-^ + A) 



aj3 



e-1 



m 



Substitute Ej and X^^ in Eq. (lASp to obtain P(t'a — )■ 1^/3). 



Appendix B: The disappearance oscillation probability of at high-energy in the 
standard three-flavor scenario in matter 



In this appendix we derive the high-energy behavior of the the oscillation probability 
P(z/^ — )■ u^) in the standard three-flavor scenario in matter. At high-energy we can safely 
ignore the contribution from the solar neutrino mass-squared difference Am^]^. In this case 



the three eigenvalues can be easily obtained as Ei, E2, i?3 = A_, 0, A+ [72[, where 

A± = ^E^LtA ± l,y(Ai?3i cos2^^i3 - A)^ + (AE31 sin 2^13)'. 

They have the following behavior at high-energy (|A£'3i| = |Am|^|/2£' ^ A): 

E, ^ AEsicj, 
E2 







E3-A, 



where cjk = cos 9jk. From these and the fact Y2'^ 
the coefficients 



AE3i|f/^3|', yi"' = AE|i|?7^3|' we get 



MM 



-13 



xr ^ 



1 



ITT 12 

^3 

Cl3 



AE^ 



31 



A 



\TT |2«2 



where Sjk = sinOjk- Hence we obtain the following behavior of the disappearance probability 
at high-energy: 



hm ' - ^('^^ ^ 



lim 



1 



E^oo {AE^i/Ay E^oc {AE^i/Af 



1 



sm 



2 cl,AE;iL 



sin^ 26. 



cf^AL 



c 

+ 4 
, 2 



13 



'13 



AE. 



31 



A 



|2„2 



2fAL\ 



23 



S23 sin^ 26*13 sin^ f 



m3| Si^sm J 
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or in other words, we have 

1 - Piu^ ^ u^) oc ^. (Bl) 

Eq. is the prediction for the high-energy behavior of the disappearance probabihty from 
the standard three-flavor oscillations, and it is perfectly consistent with the atmospheric 
neutrino data [t^. 



Appendix C: The disappearance oscillation probability of i/^ at high-energy with 
new physics: (i) and e^^ 

In this appendix we consider the disappearance oscillation probability with new physics 
in propagation given by Eq. ([2]), particularly with nonvanishing eg^^ and e^,-. To simplify 
calculations, let us normalize U£U~^ + .4. by A{1 + eee): 



USU-^ +A = A' 
where we have definecill 



[/diag 



AEoi AE- 



31 



A' 



A' 



A' 



u- 



efiJ 
I \ 



MM 



— 



A{1 
l + ee 



(CI) 

(C2) 



le' 1 


le' 1 




l^efi 1 5 







are relatively small compared with one, while \e'^. 



and |e!^^| could be of order 0{1). 

First, we evaluate the energy eigenvalues for \AEsi/A'\ ^ 1, |e'^^| ^ 1 (a = e,fi,T), 
\e'^T- — jCg^Pl <^ 1, while \e'^^\ and \e'g^\ are not necessarily expected to be small. Furthermore, 
we will ignore small corrections due to AE21/A' and sin ^13. Then the mass matrix becomes 
1 




USU 



-1 



A 



-e/i 



fie J 



"MM 



+ ^23^ 



er 

£23-523^ 



)* + C23S23?? 4r + C23^ 



e' 












e" 


e" 


MM 


M^ 








e" 



(C3) 



where 

= AE31M' (C4) 
has been introduced and the difference between e'^a and e'^a is the term proportional to rj. 



Expanding the eigenvalue equation to first order in e'^. 



and to second order in 77, 



-e/i5 fj.fj, 



e' „ and e' , we have 



-flTl 



{U£U-^ +A)-tl 
+Hl4r + K.P + |e' 



2 2 2\ 
C23"^23'7 J 



le" r + e' |e' | 

I /xrl ' tt\ efi\ 



- 2Re{e' (e;j*} + .^3(|e;j2 - e'^^)r^ - clAd 



(C5) 



^ Here, we have assumed that |1 + eee| is not very small. This assmiiption is reasonable, because otherwise 
the solar neutrino would not feel the matter effect very much and it would contradict with the solar 
neutrino data. 
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In the limit where all small parameters r], e'^^ — ICg^P, e'g^, e|^^, e^^, are zero, Eq. (ICSP gives 
the roots t = 0, 0, 1 + e!^^. If we include the corrections to first order in the small parameters, 
then we get 



t2 



l + t' 2 



± 



1 



2(1 + 



4(1 + {le^^r + 4.141' - '^^^^'eAMrT) 

1/2 



+4(141' - 4)^ - 44^'} + (1 + 44)'^ 

t3 = l + e;, + 0(r7,e;^). 
Next, we evaluate the quantities Y-^^ . They are given by 



^3 



1 



e' |2 + |e" |2 + |2 



where e'^^ are defined in Eq. (EH]). Plugging these results into Eq. ( lAllI) . we obtain 

1 



At2i(l + e;. 
-1 
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^j-^ ^1 „2 I ^/ ^ I |2 I \ji |2 I 1^// |2 



^^23^)(1 



e + e r + e r + e , 
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1 + So,e 



23tTT / / 



4 + 4^)^+141+141+1^ 



e" 1 2 

=/xr I 



(C6) 



(C7) 



where At2i = h — h and (j = 1, 2, 3) are given by Eq. flC6l) . 

Let us now consider the situation in which the limit E 00 or rj = A'm'^i/{2A'E) — 
is taken while e'^^ 7^ and e^^ 7^ are kept. In this case, from Eq. ([3]), which implies 
\e'^^\'^ < 0.1 and |e|^^p < 0.1, and from Eqs. we get 



1 - P{u^ 



At2iA'V 



4X3^^X1 -^s"") sin 



2, (l + e;jA'L 



(C8) 



where 



4^XPPXPP _ ^ 
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2Re(e^^e;,(e;,, 



1 + e' 



' ' ' J4P + 4.l4P-2Re(e' e' (e'J*) 



le' r + e' 

^1 e^il l"^^t^j 



(1 + 4. 



(1 + 4. 



At2i ~ 2t2 ^ 2 



l4l' + 4rl4l'-2Re(46;,(6^j 

1 + eL 



Recovering the original notations, we obtain Eq. ([8]) from Eq. fICSp . Since the argument 
(1 + e'^^)A'L/2 = (1 + + e^T-)y4L/2 of the sine function of the second term on the right 
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hand side of Eq. ([8]) is of the same order as in order for Eq. ([8]) to be consistent with 

the atmospheric neutrino data at high-energy (cf. Eq. JBip ). at least either of the followings 
has to be satisfied: (i) Atai = X^^ = or (ii) AX'^^X!^^ = X^^ = 0. Let us introduce the 
two quantities: 

p — 1/ |2 I \p' |2 I \J |2 
I e/j| l^/ifil ' \ iit\ 

Then, in the case of (i), we have F G 0, which imphes \e'^^\ ^ |e^^| ~ \e'^^\ ~ 0. In the 
case of (ii), on the other hand, we get 

f F \\, G 



G 



1 + 6;, 

From these equations and the fact that |e' P ^ 0.1 and |e' ^ ^ 0.1, it follows that 



i.e., in the original notation. 



2 I U |2 ^ ^TT\€-^^ 



\2 



l^e/il + |e^T| — 1 _|_ ' 
-L ~r tee 



has to be satisfied at least, if a consistent solution exists at all. The bounds from (i) or (ii) 
are stronger that of Eq. (jS]), and it is expected that the three parameters ee^, e^^, e^r are 
negligibly small compared to ege, and e,-,-, although this speculation has to be verified by 
numerical calculations. 



Appendix D: The disappearance oscillation probability of at high-energy with 
new physics: (ii) e^r - |eerP/(l + (-ee) 

In this appendix we assume that te^ and e^r are zero and discuss the case in which eee, 
and are of order 0{1). In this case the energy eigenvalue equation (]C5|) becomes 

= -t3 + ^2^^ ^ ^/ ^ ^ ^) _ ^1^. ^ ^ 4e;jr/} + 4C'^, 

where e'^p and 77 are defined in Eqs. ( 1C2I) and ( 1C4I) . respectively, and 

C'^e' -Ie' I'. 
In the case where \ri\ <C 1, the three roots are given by 



ti ~ 
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and Xf'^ are given by 



1 + f' 

1-1 ^^^23 

-T^/i/i '^3 2 _ '"3 2 



xA^t/x ^2 2 _ C 2 



The disappearance oscillation probability has the following behavior: 



As 



2 

23 



32 



C [it3 ) tl \ 2 ) 
where At32 = ^3 — ^2- By recovering the original notations, we obtain Eq. ([9]) 



Appendix E: The disappearance oscillation probability of v^^ at high-energy with 
new physics: (iii) fee? ^er and e,-,- = [ger 

In this appendix we assume C' = Crr ~ kerP = derive the high-energy behavior 
(ITT]) . In this case, introducing the new variable tp = |eer|/(l + ^ee) (cf- Eq. (fTOj) ). we have 



1 + ee 



' _ \ J |2 _ .2 
^TT — l^erl — ''/3- 



Recovering the small corrections due to ^13, we have the energy eigenvalue equation: 

= -t^ + t\l + tl + 7])-tr](l+tl - \Ue3 + Urse'j' 



•^13 "T - <-'/ -r ^13 

= -t[e-{l/cl + r^)t + r^{c'[,f/cl}, 

where c'{.^ = cos 6'{^ and 9'(^ is defined by Eq. (]F6P The three roots in this case are ti = 0, 
^2 — ^{(^i^Y 1 — l/c^- Then we have the following high-energy behavior (ITT]) in the limit 
E 00: 

1 - P{u^ ^ z/^) _ 1 - P{u^, v^) 



(AE3lM)2 (l + eee)2r/2 



■^23 


1- 


■^23 


(C'1'3)^ 1 


(C'1'3)^ 


■^23 


sin^ 


2^r3 (- 





2 



I 2 J 



sm 
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Appendix F: The analytic expression of P{i^^ — )• fe) with new physics in propagation 



In this appendix we discuss the analytic form of the osciUation probabihty Pii^fi — > t'e)- 
From the formula flA5D we have 



-4Re (X^'xr) sin 
-4Re (X'^'xr) sin 
-8Im (xfXf *) sin 



4Re [X^'xr) sin' 



AE32L 



\ 2 ^ 

f AE2iL\ . ( AE^2L\ . f AE31V 



sm 



sm 



(Fl) 



From Eq. (]A11|) Xj"^ can be expressed as 

1 



fie 



xr 



A£'2i Ai?3i 
-1 

A£'2i ^£'32 
1 



[rr - {E2 + Es)Yr''} 



(^3+^i)Fr} 
(Ei+E2)Fr} 



A^3lA^32 

where Yj^'^ (j =2, 3) is defined by Eq. flAlOp . and are given by 
¥2^^ = AE^iX^^ + AEsiXf 

Vr = {AEs.fXt + (AE^O^Xf + AAE,, {(1 + e,e)X3^^ + e.,X3^^} 

+AAE21 {(1 + eee)Xf + e.eXD . 

Here we have introduced the same quantity in vacuum: 



(F2) 



X- 



Tl Jl* 



As shown in Ref. |67|], in the limit Am^^ — )■ 0, Ej can be expressed as roots of a quadratic 
equation. First, let us review how to obtain them. With the ansatz (fT6!l . we have 



where /9 is defined by Eq. ([10 



A"^A(l + eee)/cJ, 

7 = ^arg(eer), 

and we have introduced notations for 3x3 hermitian matrices: 

f -i 




(F3) 
(F4) 



A, 



Ao 




A. = 



Av 



\ i 






(1 







j, Ag^ 








-I) 




lo 
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Here, A2, A5 and A7 are the standard Gell-Mann matrices whereas Aq and Ag are the notations 
which are defined only in this paper. 
The mass matrix can be written as 



USU~^ + A 



U£U~'e 



+ diag (A", 0, 0) 



Here, we introduce the two unitary matrices: 

U' = e^^^^e-^^^^ U 

= diag(l, 1, e'^^'^^'r^) U" dmg{e''''^^'^^ , e'""'^^'^^ , 1), 

where U is the 3x3 MNS matrix in the standard parametrization [l| and U" was defined in 
the second fine in such a way that the elements t/"^, U'J2, U'^^ be real to be consistent with 
the standard parametrization in Ref. [l[ [1. Then we have 



U£U-^ +A = e*^^'-'e-^^^Miag(l, 1, 



jargC/'., 



U"gU"-i + diag (A", 0, 0) 



xdiagd, 1, e-*^'^^^-^) gi/3A5g-i7A9_ 



(F5) 



Before proceeding further, let us obtain the expression for the three mixing angles 6'-^ and 
the Dirac phase 5" in U" . Since 



U' 



u, 



111 



u, 



/x2 



u, 



/x3 



where cp = cos (3, sp = sin/3, we get 



sin-' \U':^\ = sin-' \cpe-'^Ue3 + spe'^rs 



9';, = tan~\u:,/U:,) = tan-i {\cpe-'^U,2 + spe'^Ur2\/\c^e-'^U,, + spe'^Uri 
= ts.n-\U';jU':,) = tan-i [u^s/lc^e-'^rs - s^e'^U^^ 

6" 



(F6) 
(F7) 

(F8) 



-arg 



{U',SrKlUe2K 



As shown in Ref. [72], in the limit Amgi — )■ 0, the matrix on the right hand side of Eq. (1F5 
can be diagonalized as follows: 

U"SU"-^ + diag {A", 0, 0) - ^il 



r5»e^^"3^T^/e*^i2^Miag (0, 0, AE31) e-'^'^^^^rs"e~'^'^^^'rj,}e-'^'^^^' + diag [A", 0, 0) 



e<3^5(jiag (0, 0, AE31) e-'''^^^' + diag {A", 0, 0) 



e*^23^^r5„e^^"3^Miag (A_, 0, A+) e-'^'i^^'rj,}e-'^'^^^' 

iS"^ 



where Ts" = diag(l,l,e "^"), AE-^i = Am1^/2E, we have used the standard parametriza- 
tion [1] U" = e*^23'^7p^„g«fi3A5p-ig«6'i2A2^ g^j-^^ ^j^g elgeuvalucs A± and the effective mixing angle 

9".^ are defined by 

l(AE3i + A")±[ 
AE31 sin 2^7, 



A± = - (AE31 + A") ± -V(A^3i cos2ei'3 - A"Y + (AE31 sin2ei3 



tan2^'/3 



'13 

AE31 COS 2^i'3 - A 



The element U^2 has to be also real, but it is already satisfied because [/"j = Ut2- 
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In the present case, since we consider the oscillation probability -P(t'^ — ?■ t'e) in Korea 
at low energy, i.e., at L = 1050 km with E = 0.7 GeV, the mass-squared difference due 
to the solar neutrino oscillation gives a non-negligible contribution AE21L ~ 0.3 and the 
correction in | Ampi | / | Am ^i | becomes important. 

The discussions in Ref. 67| can be generalized to the case with non- vanishing Am^i, since 
discussions up to Eq. (IF5p are valid for a generic value of Aml^ ^ and all we have to do is 
to obtain the correction to the energy eigenvalues due to small Amg^. The energy eigenvalue 
Ej to first order in |Am2i|/|Am|i| can be computed and are given by 



^1 = A- + {A+ - (c'/a)^^"} Ai?2i 



A+ - A 
E2 = {^[^fAE^i 

^3 = A+ - -^^^ {^- - {c'[,fA'] AE,,, (F9) 

where c"^ = cos 9"/^ and sj^ = sin^^"^. As seen from these expressions, in the limit Amgi — )■ 0, 
all the quantities depend on the phase only through the combination 5 + arg(eeT-), since they 
depend only on ^'/g which is a function of the combination 5+arg(eeT)- On the other hand, the 
oscillation probability -P(f/i — )■ i^e) in Korea has the moderate dependence on 6 and arg(eer) 
separately, because the first order corrections in lAmlil/IAmg^l to Ej have dependence on 
9I2, which is approximately a function of arg(eer) only. 



Appendix G: The T violating term in Pii^fj, — > Ve) 

In this appendix, to see the contribution of the two phases S and arg(eer) in -P(t'^ — ?• i^e), 
we will study the T violating term, which is the last line in Eq. (IFip . This term contains 
the modified Jarlskog factor Im(Xf^X2^*), and it can be rewritten as 



Im(XrXf*) 



-1 



Im 



{AE2i)^AEsiAEs2 



{Yf - {E2 + Es)Y2^'} {vr - {E, + 



(Gl) 



A£'2i AE'si A£'32 

In the present case with the ansatz (fT6l) . instead of using the explicit expressions f lF2p . it is 
convenient to work with the following form: 

{USU-' + A),e = {USU-'),, 



-'2 
^3 



[U£U-^ 



A 



{US'U-')^e + {USU-')^eAee + {U SU 



r-1' 



A 



SO that the factor \m.{Y'^^ {¥2^)*} in Eq. (IGip can be expressed as 

Im{F3'^^(F2'^^)*} = Im [{U£U-Xe {{UE^U"'),, + {U£U-'),rAre] 

= \rn{Yr{Y2'r}s., + im{r3''^(rr)*}NP , (G2) 

where the term {U£U~^)^eAee in Y^^ has dropped as in the case of the standard scheme, 
because \{U £U~^) ^e\'^ Aee is real. In Eq. ( ]G2p we have introduced the notations: 



Im{Kr^(F2'^^)*},,, ^ Im \{USU-Xe{U8'U-'- 



lie 
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Im 



= AE2iAE,iAE,2lm{Xt{Xn*} 
is the Jarlskog factor in the standard three- flavor scheme and 



(G3) 



Im{r3^'=(Fr)*}NP = lm{{USU-XeiUSU-'),rAre} 

= Im [AieerTiAEsiiXn* + AE,,{Xn*}i^E,,X^^ + AE,,Xn] 



{AE,,flm{{Xn*Xfe''''} 



(G4) 



is the extra contribution to the Jarlskog factor due to new physics. If ^13 is small, then the 
dominant contribution in the new physics term (lG4p comes from the middle one which is 
proportional to yl|eeT|Ai?3iAii^2i|-^2'^||-^3^|, and this should be compared with the standard 
factor ( 1G3I) to examine which contribution dominates the T violating term. 
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